High-order perturbations of a spherical collapsing star 
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was developed. In this article, we apply it to the particular case of a perfect fluid background. 
We have expressed the perturbations of the energy-momentum tensor at any order in terms of the 
perturbed fluid's pressure, density and velocity. In general, these expressions are not linear and have 
sources depending on lower order perturbations. For the second-order case we make the explicit 
decomposition of these sources in tensor spherical harmonics. Then, a general procedure is given 
to evolve the perturbative equations of motions of the perfect fluid for any value of the harmonic 
label. Finally, with the problem of a spherical collapsing star in mind, we discuss the high-order 
perturbative matching conditions across a timelike surface, in particular the surface separating the 
perfect fluid interior from the exterior vacuum. 
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I. INTRODUCTION 



The detection of gravitational waves is currently con- 
sidered to be one of the most important open prob- 
lems in astrophysics and experimental physics. Aside 
from providing a test of the theory of general rela- 
tivity in the strong-field regime, it will open a new 
window for astrophysical observations, giving rise to 
the era of gravitational-wave astronomy. To this end 
the ground-based laser-interferometric detectors LIGO, 
VIRGO, GEO600 and TAMA 0-Q are now collecting 
data at or near design sensitivity, while next generation 
ground-based detectors like the Einstein Telescope (ET) 
are in the research and development phase [7] . The space 
based interferometer LISA [8J, scheduled for launch in 
the next decade, will complement such observations with 
high precision in the low-frequency range. Gravitational 
wave detection and interpretation requires a close inter- 
play between the experimental side and theoretical pre- 
dictions as physical signals need to be dug out from a 
noisy data stream. Various astrophysical scenarios are 
expected to generate gravitational radiation that will be 
detectable on Earth. Most prominently among these fea- 
ture the inspiral and coalescence of compact binaries, 
black holes and/or neutron stars, as well as stellar col- 
lapse. The work of this article is primarily addressed at 



the modelling of the latter type of sources although we 
believe our results to be of interest for a variety of ques- 
tions pertaining to gravitational wave source modeling. 

Simulating neutron stars and black holes in the frame- 
work of fully non-linear general relativity is only possible 
using numerical methods on super computers, an area 
of research commonly referred to as numerical relativity. 
This field has made enormous progress recently; see, for 
example, the reviews [9l-fllj]. In spite of this progress, 
approximative techniques such as the post-Newtonian 
expansion and perturbation theory still play a crucial 
role in the modeling of astrophysical sources of gravita- 
tional waves. First, numerical simulations are compu- 
tationally too expensive to cover more than a few tens 
of orbits in inspirals, so that complete waveforms are 
now constructed by hybridization of numerical and post- 
Newtonian results [TH [l3[ or calibration of the effective- 
one-body model via numerical relativity simulations [l4l - 
Ojj]. Second, other types of sources completely elude a 
full numerical treatment. For example, simulations of the 
so-called extreme mass ratio binaries are practical only 
in a perturbative framework; see fl7l. [l8l a nd references 
therein for a discussion as well as for the most 

extreme mass ratios simulated in numerical relativity. 
Perturbation theory also provides a natural framework 
for interpreting numerical results, as for example in the 
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case of black- hole and neutron-star oscillations [23|, |24j . 
Finally perturbative calculations efficiently provide re- 
sults with high precision even when involving numerical 
methods and are thus ideally suited for the study of a va- 
riety of astrophysical and cosmological problems as well 
as benchmarking the accuracy of full blown numerical 
simulations. 

Purpose of this work is to provide a general gauge in- 
variant framework for perturbation theory of arbitrary 
order on time-dependent and spherically-symmetric per- 
fect fluid backgrounds. This work constitutes apartic- 
ularization of the formalism developed in Refs. [J by 
including matter obeying perfect fluid equations of state. 
In particular, we will also include a careful analysis of 
the matching conditions necessary to obtain a consistent 
treatment of the stellar interior and the vacuum exterior 
spacetime. 

Static backgrounds represent a particularly simple sub- 
class of configurations because they allow for a Fourier 
expansion of the perturbations in such a way that reduces 
the problem to a set of ordinary differential equations de- 
pending on radius only. Examples for the general type 
of background configurations considered in this work in- 
clude radially collapsing or exploding stars. The collapse 
of homogeneous dust spheres studied extensively in a se- 
ries of papers by Cunningham, Price and Moncrief [25l — 
H3] can be viewed as the special case of a perfect fluid 
with vanishing pressure. In this case the different layers 
of the star do not interact which leads to a simpler type 
of motion but may result in unphysical consequences. 

The perturbations can be either radial or non-radial, 
although in some scenarios it turns out to be convenient 
to absorb radial perturbations in the spherically symmet- 
ric background configuration. See, for example, Ref. [28[ 
for a fully non-linear numerical treatment of radial pulsa- 
tions of neutron stars. Non-linearities are naturally incor- 
porated in perturbation theory by including higher-order 
terms in the calculation. First-order perturbation theory, 
for example, linearizes the Einstein equations in the first- 
order metric perturbations whereas second-order 
theory results in a set of differential equations linear in 
the second-order metric perturbations <2} ft. MI / but con- 
taining source terms quadratic in the ll} h^ u . This hi- 
erarchical structure can be extended to arbitrary order 
although the equations quickly become overwhelmingly 
complicated. 

The application of perturbation theory to compact 
stars in a general relativistic framework was pioneered 
in the 1960s by Chandrasekhar's analysis of the stability 
of static-background models against radial perturbations 
[2!j[3(|- Next, Thorne and collaborators established the 
theoretical basis for the stud y o f non-radial perturba- 
tions of perfect fluid stars [3lH35j . Further development 
and applications can be found in Refs. [36M39j | . Time- 
dependent background configurations were first studied 
in the above mentioned work by Cuningham et al. I2E 
on collapsing dust spheres. Seidel and coworkers f.4C 
applied the formalism of Gerlach and Sengupta [4c 



to generalize the work of Cunningham et al. to gen- 
eral perfect fluid equations of state. Motivated by this 



work, Gundlach and Martin-Garcia [45[ developed a co- 
variant and gauge-invariant framework to analyze arbi- 
trary first-order perturbations of a spherical perfect fluid. 
This framework was later used by Harada et al. [46| to 
analyze axial perturbations of stellar collapse. 

To our knowledge higher-order perturbation theory of 
stellar models has so far only been applied to the mod- 
eling of slowly, uniformly rotating stars. The key idea 
here is to incorporate the rotation in the form of a first- 
order axial perturbation. This idea has again been pio- 
neered in the Cunningham et al. series 27J in the spe- 
cial case of a collapsing dust sphere. Applications to a 
wider class of perfect fluids have so far been restricted to 
static background models [13, |48| . Differential rotation 
has been modelled in the context of the Cowling approx- 
imation 1 [49|, [5(| ■ More recently it has become possible 
to study perturbations of fast rotating axisymmetric con- 
figurations, still in the Cowling approximation [5lTj53j . 

Our work represents the natural extension of the ex- 
isting literature by constructing a general framework for 
second-order perturbations on a time-dependent, spher- 
ically symmetric background spacetime containing mat- 
ter in the form of a perfect fluid star. The starting point 
for this construction is given by the second-order gen- 
eralization of the Gerlach-Sengupta formalism developed 
in Refs. [l|, @|. Specifically, we will apply this gauge- 
invariant formalism to the case of a perfect fluid star fol- 
lowing the notation and techniques presented in Ref. [45j . 

A key ingredient in our study is the matching at second 
perturbative order of the stellar interior and the vacuum 
exterior parts of the spacetime. Its first order analogue 
was analyzed in a covariant and gauge invariant manner 
in [54[ , but we will see that the extension to higher orders 
represents a non-trivial problem. Its solution requires a 
careful analysis as to which quantities must be continuous 
across the stellar surface. Our analysis of this particu- 
lar point will be valid at any perturbative order and we 
believe, it will also be helpful in clarifying the matching 
conditions in certain background configurations. 

A brief overview of the quantities and the equations 
governing their behavior is given as follows. The back- 
ground metric in the stellar interior is described by the 
scalars U and W defined in Eq. (f2"5]l and the structure 
coefficients fi and v given in Eq. (|2"Tj) , all defined with 
respect to the background velocity field. The fluid is de- 
scribed by its density p, the entropy per fluid element s 
and the pressure p. These variables are determined by 
Eqs. (|2"4")) - (pri"j) and the equation of state that provides p 
as a function of p and s. 

The perturbations at order n are naturally divided into 
axial and polar modes and further decomposed into mul- 




1 In this approximation one neglects metric perturbations and thus 
evolves fluid perturbations only. 



3 



tipoles with labels I and m. For simplicity, we will ommit 
these labels from the perturbation variables in this sum- 
mary as is also done below in the derivation. The axial 
perturbations are given by the velocity component /3 (|41[) 
and the Gerlach-Sengupta master scalar II (|68[) which 
are evolved according to Eqs. (|73|) and (|64|) . The metric 
components <5 and A (|61[) can be reconstructed from the 
master function II according to Eqs. (T7UI) . (fTTj) . The po- 
lar perturbations are given by the velocity components a 
and 7 (|41l) the energy and entropy perturbations u) and 
a (|3~4")l . (|3"5j) . the metric perturbations rj, x an d "0 de- 
fined in Eqs. $2$, (JB3J and the scalar £ (|I2J|. This set of 
perturbations obeys the evolution system (|9Tj) , (|92|) - (|97| . 
(HQ6]). 

The remainder of this article is organized as follows. 
Section [TT] presents a brief summary of the generalized 
formalism introduced in Ref. 0,0. In Sec. EH 

we in- 
troduce the notation for the background fluid variables 
and their perturbations. We further express the pertur- 
bation of the energy-momentum tensor in terms of these 
variables. The evolution equations resulting from the 
perturbed Einstein and matter equations are presented 
in Sec. IIVI The high-order matching conditions are the 
subject of Sec. [V] This section is quite independent from 
the previous ones, so the interested reader is referred di- 
rectly there. We conclude in Sec. I VII with a discussion of 
future applications and extensions of our work. 

The intensive tensor computations in this work have 
been performed with the Tensor Computer Algebra 
framework xAct (55|, based on Mathematica. Of particu- 
lar importance has been the package xPert [56[ for high- 
order metric perturbation theory around curved back- 
grounds. 



II. HIGH-ORDER GERLACH AND SENGUPTA 
FORMALISM 



where gAB is the metric on A4 2 , ^y a b is the standard met- 
ric on the unit sphere and r, Q are scalar fields on M. 2 . 
We denote covariant derivatives associated with the dif- 
ferent metrics by 

9^u-\ = 0, 9ab\c = 0, lab:c = 0, (3) 

and define the vector va = t\aI t - The Einstein equations 
for the metric ([TJ with matter tensor ([2]) are given in [43[ . 

Tensor fields of any rank s on the sphere will be decom- 
posed using a basis of tensor spherical harmonics labeled 
by multipole indices / and m. Such basis can be con- 
structed from the symmetric trace-free (STF) tensors 

Zr ai ... as = (!T:ax...aJ STF , (4) 

•y-m b ym fr\ 

■A/ ai...a s — e (ai 6a 2 ...a s )j l°J 

which are respectively polar and axial, together with the 
metric ^ a b and the antisymmetric tensor e a b For the 
particular case s — 0, those objects must be read as Zf 1 = 
Y™ and X™ = 0. 

At first order modes with different I, m and polarities 
decouple. However, at second and higher perturbative or- 
ders, the non-linear coupling between first-order modes 
results in products of tensor spherical harmonics with dif- 
ferent labels (l,m,s) and (/,m, s). Those products can 
be decomposed into a linear combination of harmonics 
(/, m + m, s + s) with an explicit formula involving coef- 
ficients 

where C stands for Clebsch-Gordan coefficients. These 
^-coefficients encode the geometric selection rules that 
determine which pairs of modes do actually couple. See 
1] for full details. 

B. Non-spherical perturbations 



A. Background spherical spacetime 

The spherically symmetric nature of the background 
spacetime enables us to naturally decompose it as A4 = 
M 2 x S 2 , where M 2 is a two-dimensional Lorentzian 
manifold and S 2 the unit two-sphere. In order to distin- 
guish tensor fields residing on these different manifolds, 
we will use Greek letters (/x, v, . . .) for four-dimensional 
indices, capital Latin letters (A, B, . . .) for indices on M 2 
and lower-case Latin letters (a, b, . . .) for indices on the 
sphere. With this notation, we can decompose the back- 
ground metric and energy-momentum tensor as 

9nv{x D , x d )dx fl dx 1 ' — gAB{x D )dx A dx B 

+ r 2 (x D ) lab {x d )dx a dx b , (1) 
t^u(x D , x^dx^dx 11 = tAB{x D )dx dx B 

+ \r 2 (x D )Q{x D ) lab (x d )dx a dx\(2) 



Perturbative calculations are often formulated in terms 
of a one-parameter family of spacetimes (Ai(e), g(e)). 
The expansion parameter represents a measure for how 
strongly a physical system deviates from the background 
configuration corresponding to e — 0. In complete anal- 
ogy, any tensorial quantity of the system is represented 
by a family f2(e), and expanded around its background 
value Q = 0(e = 0) according to 

^) = "+E £ "^P- (?) 

71=1 

A key feature of this expansion is that the perturbations 
A n [f2], also denoted by { ™ } f2 in this work, are tensors on 
the background manifold. We have thus converted a sin- 
gle problem on an unknown manifold into an infinite hi- 
erarchy of problems on the chosen background manifold. 
There remains, however, one important problem; we need 
to specify a diffeomorphism 4>(e) relating points on the 



4 



background A4 and the perturbed manifolds A4(e). The 
actual values of the perturbations depend on this 

choice of gauge. 

We next explicitly decompose tentorial perturbations 



in terms of the basis of tensor harmonics introduced in 
Eqs. (UJ) and ([5]). In this expansion the metric perturba- 
tions are given by, 



A n [g^] ee <»>ft. 



E 



I AB 



I A^i b 



'"'ff(™B Zp a + ^h? B Xp a ^K™ r 2 lab Z™ + <">G[ 



^hf A X™ b 
■ r 2 Z™ ab + 



X 



l ab 



(8) 



At each perturbative order n, the gauge freedom enables 
us to choose freely four of these functions, three polar and 
one axial. Alternatively to thus fixing the gauge, we can 
construct gauge-invariant combinations of the perturba- 
tion functions and work with those quantities. We will 
follow this latter approach, though with a construction 
based on the choice of a particular gauge. 

For this purpose we consider the most natural gauge 
choice in spherical symmetry: Regge- Wheeler (RW) 
gauge [57j |. It is defined by setting 

{n} Hf A =Q, M Gr=0, { "V = 0, (9) 

for all n > 1 and I > 2. We note that this does not con- 
stitute a rigid choice, that is, there remains some gauge 
freedom in the I = and I = 1 modes. Next we use 
an important result from Ref. [2j, namely a procedure 
to define gauge-invariant quantities whose values coin- 
cide with those of the perturbations in RW gauge. The 
key idea of this procedure is to employ the formula for 
general gauge transformations parametrized by the gen- 
erators changing the perturbation A™ to 

n . ^ 

A n [Ql + V - V - 

[ J ^( n - m )!^ 2! fe 2...(m!) fe ™/fc 1 !...fc m ! 

m=l v ' (K m ) 
I 



x c% r x\z k A n - m [n}, (10) 

where the second sum extends to the m-tuples 

{K m ) = Uki,...,k m )€nZ-, X ./A'.. '"j- (11) 

N being the set of non-negative integers. We then need 
to choose the vectors <m} £ such that they take the metric 
perturbation (jSJ) from an arbitrary gauge to the RW form 
©. Naturally, the vectors {m} £ depend on the perturba- 
tions m hn V themselves, and so the expressions for the 
gauge invariants are highly nontrivial beyond first order. 
Explicit formulas at second order, already expanded in 
terms of the harmonic coefficients, are given in [2j. 

Given that the perturbations in RW gauge and the as- 
sociated gauge invariants coincide in value, we will per- 
form all our computations in RW gauge, for simplicity, 
and then use the expressions in Q whenever we want to 
generalize the results to an arbitrary gauge. 

In this way, the decomposition of the perturbations of 
the metric and the energy-momentum tensor will be as 
follows, 



l AB ^1 
m 
I a 



{n} m 



^Kf B X 



K?A Xp b 

M /Cr r 2 lab Z? 



(12) 



A*[t 



E 

1. rn 



l B &i a 



I AB 



i^n A Z rn b + X 



lab' 



J l ab 



rn 
I b 



X! 



(13) 



Note that polar (axial) harmonic coefficients are repre- 
sented by upper-case (lower-case) letters. In order to 
agree with the construction by Gerlach and Sengupta 
[ij, |44|, we combine the invariants { " } V>[™4 an d { ™ >K z™4 
to form a new matter invariant quantity 

{ "Vi"A-§ { " } «rA. (14) 



III. PERFECT FLUID 

The decomposition of metric and matter perturbations 
given above is valid for arbitrary types of matter sources. 
In this section we discuss in detail the special case of 
a perfect fluid. We first formulate the equations for 
the background spacetime in Sec. IIII A[ following [45| . 
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Then we address the description of fluid perturbations in 
Sec. IIII Bl generalizing the notations and results in [45[ 
from first to higher orders. The metric and fluid pertur- 
bative evolution equations will be analyzed in the next 
section. 



A. Background perfect fluid 

The energy-momentum tensor for a perfect fluid with 
four velocity w M , total enery density p and pressure p is 

tnu = (fi + p)u^i v + pg^u- (15) 

This corresponds to the special case Q = 2p in the general 
expression ((2]). We will use a generic equation of state 
for a perfect fluid, in the form p = p(p, s) where s is the 
entropy per particle. It is convenient to introduce the 
partial derivatives of the pressure as additional quantities 
according to 

Here, c 2 is the adiabatic speed of sound, but we are not 
aware of a comparatively intuitive meaning of the second 
quantity. 

In spherical symmetry the four-velocity of the fluid 
takes the form = (u A , 0). This vector defines a unique 
outward pointing spacelike unit vector ua = —€abu b on 
ftA 2 , with €ab being the antisymmetric volume form on 
A4 2 . These two vectors form an ortho normal basis on 
this manifold, which can be used to decompose all geo- 
metrical objects, for example 

9AB = —UAUB + nATlB, £AB — flAUB ~ UATT-B- (17) 

Using this decomposition, the M 2 part of the energy- 
momentum tensor @ takes the form 

t A B = PUaub +pn A n B - (18) 

For scalar fields / on M 2 we further define the frame 
derivatives 

f = u A f\ A , f = n A f lA - (19) 

A straightforward calculation shows that the frame 
derivatives obey the commutation relation 

(/)'-(/')'=M/'-^/, (20) 

where 

fi = u A {A , v = n A \ A - (21) 

These "structure" functions are the components of the 
covariant derivatives of the frame vectors in their own 
frame 

ua\b = n A (n B p-u B v), n A \B = u A {n B p-u B v). (22) 



We are now in the position to formulate the Einstein 
equations exclusively in terms of the background scalars 
p, p, p, v and 

U = u A v A , W = n A v A - (23) 

For our spherically symmetric background we thus ob- 
tain expressions with remarkable symmetry for the four 
independent field equations 

U' =W{p-U), (24) 
W =U{v-W), (25) 

W' = -471-p- W 2 + Up+ -g-, (26) 

U = -4irp- U 2 + Wv~ -3, (27) 

with the Hawking mass [58| in spherical symmetry 
(Misner-Sharp mass) 

M = T - (1 - r. A r A ) = V - [1 + r 2 (U 2 - W 2 )} . (28) 

We still need to consider the equations of motion for 
the perfect fluid. Conservation of energy-momentum re- 
sults in 

P +(p + p)(2U + p) = 0, (29) 
c 2 s p' + pCs' +{p + p)v = 0, (30) 

which represent energy conservation and the Euler equa- 
tion, respectively. Finally, a perfect fluid does not dissi- 
pate energy and hence the entropy of each fluid element 
is conserved, 

s = 0. (31) 

The system (|24ti31|) fully describes the dynamics of our 
background spacetime. In our derivation of the perturba- 
tion equations below, they will further enable us to sim- 
plify those coefficients which depend on the background. 

B. High-order perfect fluid perturbations 

Our formulation of the fluid perturbations proceeds in 
two stages. First, we expand the perturbations of the 
energy-momentum tensor in terms of perturbations lni p, 
{ " } s and of density, entropy and four-velocity, re- 

spectively. Second, we need to expand the latter as series 
in appropriate tensor spherical harmonics. 

The explicit form of the n th perturbation of the energy- 
momentum tensor of the perfect fluid (|15l) is obtained by 
using standard combinatorial formulas [56j. The result- 
ing expression is given by 

= A" m = E ( I ) { w p { " v 

+ ll{ n ] Z k k )( ik} P+ W P) VV}: (32) 

j=k 
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with {n) u^ = A n [Ufj]. This can be further expanded 
along the same lines, using a high-order generalization of 
the chain rule on the equation of state, 



{«>, 



2! fe 2...n! fc "2! r 2...n! r «fci!...fc„!ri!...r„! 

where the sum is restricted to the following 2n-tuples, 

{ki,...,k n ,n,...,r n ) € N 2 n ;Y,j( k j + r j) = 

i=i 

and we have defined K = X)j=i kj and R = Y^j=i r j- ^ n 
practice derivatives of the pressure will be replaced by the 
sound speed c s and thermodynamic factor C, defined in 
Eqs. (1161) . The expression for the pressure perturbation is 
completed by the harmonic decomposition of the density 
and entropy perturbations 



i .in 



(34) 
(35) 



Our formulation of the perturbations of the four veloc- 
ity is guided by the normalization condition u^u^ = — 1, 
which holds at any perturbative order, so that 



A> p5 ^] = 0. 



(36) 



By applying the Leibniz rule n times and separating the 
terms linear in perturbations of order n, we can rewrite 
this equation as 



-2^u^ - u.u u A n [g^] 



EE 

/c— 1 i—k 
n-l 

+E 



k\{i-k)\{n-i)\ 



(37) 



On the other hand, the 71 th order perturbation of the 
inverse metric is given by 

M ft-m- 

x ( fe i>^ a ^hj... {k ^ } h a p ffcml V, (38) 

where this sum extends over all sorted partitions of n 
such that k\ + . . . + k m = n with m < n. In spherical 
symmetry the background four-velocity has no angular 
components (u a = 0), so that the constraint (|37|) only 
contains non-angular components of the n th perturbation 
of the four-velocity. Since the quantity { ™ } T defined in 



Eq. (|37[) contains only perturbations up to order n — 1, 
we can make the following ansatz, 



1 



1 



nrAB- 



A" [g 



U B 



in A , (39) 



that reproduces the first-order formula used in Ref. [45J 
and trivially satisfies the constraint (|3T[) . In this way, the 
three independent components of the perturbations of 
the four-velocity are encoded in the function { ™ } 7, which 
will be expanded in a harmonic series with coefficients de- 
noted by { ">7™, and the two angular components {rl} u a , 
with harmonic coefficients and 

The object { ™ } T, polynomial in lower-order perturba- 
tions and vanishing for n = 1, is an example of a high- 
order perturbation source. These sources, whose compu- 
tation is one of the main objectives in this article, have 
all a similar internal structure, which we now illustrate 
at the second-order perturbative level. 

We begin by decomposing the four velocity at first or- 
der (cf. [45j). 



= E 



I, m 



I ll} n n AB u B ) zr n dx A 



( {1} a r i n Zra + {li (3rX 7 i n a)dx a . (40) 



We recall our discussion below Eq. ^ and emphasize 
that all terms in this equation can be considered as gauge 
invariant. 

At second order, we obtain quadratic first-order terms 
which arise from the right-hand side of Eq. (pi!?)) and the 
perturbation is given by 

^dx^ = J2( m TnA + \ {2} 1CTabu b 



mmm 

i i i A 



LI 



Zl n dx A 



l2} a? l Zl n a + {2i (3rXl n a yx a . (41) 



Here we have defined the alternating sign 
e= (-l) r+ '- ! , 



(42) 



which is a function of the labels I, I, I, even though we 
ommit these in the notation, and the sum 



E-EE E E 

[=0 1=0 m=-l m=-i 



(43) 



with the restrictions \l — l\ < I < \l + l\ and m + ih = m 
for given I and m. The source terms {ty U can be classified 
into two groups; those of the form {+ t4 are formed by 
polar x polar or axial x axial terms, whereas the terms 
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( XI consist of mixed polar x axial terms. This structure 
becomes clear in the explicit form of the sources, 



*i i i A 



77 — A fK,Bc n u 



b„,c 



1 



K. aL, KcDU B vy -K^Kcdubu 
u B K, BC K AC 

\rhl\ U A 



aa + (3(3 + R,bu b kqu 

k B u B k A 1,(44) 



B..C 



*l i I A 



Zkbu (3 — kbkcu u 
^E^I^uaI&P + ak B u B }■ f 4.", i 



Here, the carei and bar denote first-order harmonic com- 
ponents with harmonic labels (l,rh) and (i, to), respec- 
tively. Even though these individual source terms are 
not symmetric under the interchange (l,m) <H- (l,fh), 
their sum of type (|43[) becomes symmetric in (I4ip . In 
order to save space we will follow this notation for the 
explicit form of source terms quadratic in the first-order 
perturbations throughout the remainder of this work. 

From now on we will simplify our notation for the first- 
and second-order case by dropping the labels n = 1 and 
ri = 2. The perturbative order will become clear from 
the context in the following equations. Note also that the 
expressions for the first-order perturbations can always 
be reconstructed from their second-order counterparts by 
dropping the quadratic source terms. 

With the perturbations of the four-velocity given by 
Eqs. (|40|) and (1441) . we can perform the harmonic decom- 
position of the second-order perturbation of the energy- 
momentum tensor (|32l) . The axial components are given 
by 



VT4 = (3T{p + p)u A ~i 



\)crnr\ 

c i i 



mmm 

i i A ' 



(46) 



;)pm m m 

C T i i ' 



(47) 



where E are first-order quadratic sources that will be 
given below. Similarly, we obtain the polar part 



AB 



(P + P) 
1 



jriuATiB + n A u B ) 



-m n AC u B + K? BC u A ) 



u c + P K,T AB 



puj^(u A u B + c 2 s n A n B ) + Cpa^n A n B 

E{^)cmmm 
>£[ l l AB, 



u 



{p + p)aTu A + ^ ) £rtTA, 
pK,T + c 2 s p^ + CpaT 



mmm 
III' 



U 



m rti rn 

i i i ■ 



(48) 
(49) 

(50) 
(51) 



We note that the polarity sign, that appears as a left su- 
perindex of the axial sources is — e, whereas that of the 
polar sources is e. Also, some source terms appear in 
the same form in the polar and axial equations, as for 
example on the right-hand sides of Eqs. (141)1) and (|4U)) or 
Eqs. (|4"T)) and (fB"Tj) . This is a general feature of quadratic 
sources in second-order perturbation theory; see for ex- 
ample Q. 

We conclude this discussion by giving the sources 
explicitly in terms of first-order perturbations and E- 
cocfficicnts, 



1,1 



( + )cmmm 

I l AB 



-2(p + p)E 1 ll g> l {k C K {A u B) uc - u A u B aa + (3(3 -2(3k c u c } 

EftffiA 2pc 2 s QjK AB + 2p7 [Ca +u(l + c 2 s )] {n B u A + n A u B ) + 2C 'pb(fZ A c : n B n c + JC c[B u A] u c ^ 



+ 2pu (1 + c 2 s ) t c(B u A) u c + {p + P ) 



2^(n A n B + u A u B ) - 2K. K, D ( A u B )U C 



1 „ _ 3 „ 

2jfcc(A n B)U C - 2^KcDn c u D u A u B + -JC A cu c JCbdu d - -ICc f ICdfu c u d u a ub 

( „dC n „dc 2 „dc 2 
1 2uj<j-z Jl + pww 



(-)cmm m 
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y ds ds dp 



l I l AB 



2 E l U iAl a ( kCu C ~0){P+ P)U A U B , 



(52) 
(53) 
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2iS H^(p + p) 



in A + -K-abu 



in A + -K-abu 



+ pu)(l + c s )ua + CpauA >, 



4/ 



'I i i 



2E± l t % l (aa-f3/3)(P + p), 
if 



/9 + pj(3[(£^ + l)w + C(t]ua + p{Ca + c 2 <Jj)R a 

2 - _ f -1 -_ ( Br 2 1 BC n Br 2 ~\ 

^(p + p)E^%{aa + pp) + 2pE° %^C{u>cl + aC) + ua^ + -<t<j— + P - uju>-^ j, 



IV. SECOND-ORDER EVOLUTION 
EQUATIONS 

We now turn our attention to the evolution equations. 
The results derived in the previous sections facilitate the 
formulation of the evolution equations at arbitrary or- 
der. In order to keep the complexity of the equations in 
this article at a manageable level, however, we will re- 
strict this discussion to the simplest nonlinear case, the 
evolution equations for second-order perturbations. 

Note again that our calculation represents a general- 
ization of the first-order results of Ref. [45| . Because the 
equations at second order share the linear part with those 
at first order, our main task will be the computation of 
the additional source terms quadratic in the first-order 
perturbations. 

In the Einstein equations these quadratic source terms 
contain two parts; (i) terms arising from the second-order 
perturbation of the Einstein tensor which were denoted 
by wgpp™ in their original derivation in Ref. [lj; (ii) 
sources w^ppm^ obtained from the perturbation of the 
energy-momentum tensor, which have been given explic- 
itly in Eqs. (|52H59[) for the case of perfect fluid matter. 
It is convenient to combine those two types as 

III = III III ' V DU J 

for all tensorial, vectorial and scalar sources. The per- 
turbative evolution equations for the matter can be com- 
puted straightforwardly by direct perturbation of their 
nonlinear counterparts. For simple systems, including 
perfect fluids, they can also be obtained from the equa- 
tions of energy-momentum conservation. In particular, 
the quadratic sources for the fluid equations of motion 
can be computed from the sources of energy- 

momentum conservation defined in Eqs. (100-102) of Ref. 

a- 

Once again, in order to simplify our notation, we will 
omit from now on the harmonic labels {I, m, I, m, I, fh} 
from the source names and merely write {ry P, <cy PA, 
{iY PABi and so on. We also omit the labels {I, m} from the 
second-order perturbations. They will be distinguished 
from those at first order because the latter are always 
denoted with hat or bar overscripts. 



(54) 
(55) 
(56) 

(57) 

(58) 
(59) 

I 

Let us summarize the variables that will be used to 
describe the different perturbative degrees of freedom. 
The fundamental variables describing density oj, entropy 
a and velocity perturbations {a, /3, 7}, have already been 
defined in Eqs. ([51]). ([551) and 611) respectively. All 
other tensorial perturbative variables will be decomposed 
in the natural background frame provided by the fluid 
velocity (u A ,n A ) as follows. We encode the axial metric 
information in the scalars {<5, A}, defined as 

ka = Sua + Xua, (61) 

whereas the polar tensor ICab is decomposed in three 
scalar functions, 

£ab = r)(n A n B - uau b ) + 4>(uau b + n A n B ) 

+ ^(uatib + n A u B )- (62) 

Note that the scalar ip introduced here, and that will be 
used in the rest of the article, is different from the har- 
monic function ip in the decomposition of the perturba- 
tions of the energy-momentum tensor (|13[) . In addition, 
for future convenience, the metric component function <f> 
will be replaced by the new variable x> defined by 

4>= x + IC-v- (63) 

Finally, the fourth polar geometric degree of freedom will 
be described by the scalar gauge- invariant /C. 

Regarding the linear evolution equations, we sill sepa- 
rate the polar and axial parts for arbitrary (second-order) 
harmonic label As shown in Ref. [45j . for those cases 
with I > 2, the axial sector reduces to two independent 
equations, one for the axial gravitational wave (II, de- 
fined below) and the other for a fluid perturbation (/3) 
describing rotation. On the other hand, the polar sector, 
also for I > 2, contains four fundamental equations for 
four variables, respectively describing the polar gravita- 
tional wave (x), the sound wave (K,), the entropy per- 
turbation (a) and non-rotational tangential fluid motion 
(ijj). All other metric and fluid perturbations can be re- 
constructed from those using the constraint equations. 
In particular, the function p turns out to be vanishing 
in the case I > 2. The particular cases / = 0, 1 are 
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non-radiative and the corresponding equations present a 
different structure. In particular it is not possible to con- 
struct gauge invariants for I = OA, and we will need to 
resort to adequate gauge fixing [3l|, H3] . 



A. Axial perturbations 

1. The case I > 2 

There is a single axial fluid perturbation P (defined in 
Eq. (|4Tj)). which obeys the following transport equation, 
equivalent to the axial part of the perturbed equation of 



energy- momentum conservation (cf. Eq. (99) of Ref. [l|), 



\Q* + 2U)P = i E ( ~ £) ^ 



E 



f-.W (I-l)g + 2) f _, t 



LI 



1 r 



,2(-e) 



2r 2 



(64) 



The source term ^B is given in expanded form as follows: 



<-> B . .Mll^- 

- i E^fjp'fr (^ + A (27 + 4 ) + -^-((S + <&> + C*) + 2f + 4' + * + 3« 

+ 2^-^{p + 2U)(pCa + pc 2 s uo - puo) + 2W{2^ + 2(ji + 2U)(a^ + up^)) \, (65) 

P + P OS Op J J 



(+) _ (/-!)(/ + 2) lTm 



E, 



lrnlK aa- PP) 

l\ti\ 1 ^±^2P5 + 2^j - + C °) + T £ -S 2f l~ *M1 + & - '■'( 1' > ' ) 



P+P P+P 

2 P , 2- , , 1 



a(2f + + &({x ~ 2r/ + K) + —^{u + c 2 s lo + Ca)) + — af — ^ + A (2 7 + 



+ c?a(^ + 2C/)(/C + x - 2?7) - 2at2>— ^— (1 + c?) - 2aC&—?— + 2a(p + 2U)u-^—(l + c?) 

p+p P + P ' P + P 

^ p(l + 3c.) _ _^ + _ 2 _ (2 . + ^ + pCV ( _ _ 2 _ + ^ _ 2 _ ( ^ + 2f7 ^ pCa 



P 



P + P 

2a(» + 2U)(*^+u>p^] 



P + P 



P + P 



(66) 



In order to simplify this expression we have assumed that 
the first-order harmonic numbers I and I are greater or 
equal to 2 and, thus, both f\ and 77 are vanishing, as will 
be made explicit after Eq. (|9"Tj) . We could expand the 
time derivatives of the first-order perturbations a, u> and 
a, using their evolution equations, but it does not lead 
to any simplification. 

At first order the axial perturbation P evolves freely 
on the background, in the sense that it obeys a transport 
equation in which the axial gravitational wave does not 
appear. We see in (|65|) that the source terms ^B extend 
this fact to second order because they do not contain any 
axial metric perturbations, due to a geometrical cancel- 
lation of the ^-coefficients of 6t8. Such a cancellation 
does not occur for the source terms in (|66[) . which 



contains a term 

^fe^M. (67) 

This allows a coupling between the axial gravitational 
wave 5 and a rotational mode p, generating a second- 
order rotational mode P for odd 1 + I — I. Note that in 
such a case we have SjUj?/ = —E^l^i (see [l|). 

At any given perturbative order n the axial gravita- 
tional wave is fully described by the Gerlach-Sengupta 
master scalar 

{n}n _ e AB ( __^A (68) 

V r J \B 

For I > 2 one of the perturbed Einstein equations pro- 
vides a simple way to reconstruct the vector ka from the 
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master scalar II: 

(1-1)(1 + 2)k a = 

= 167rr 2 ^ - e AB (r 4 U) lB + 2ir 2 ^ ( ~ e) S A , (69) 

1,1 

(cf. Eqs. (80) and (D5) of [J), and hence there is a one- 
to-one correspondence between II and ka, assuming the 
matter perturbations and sources are known. Expanding 
this relation in the natural frame, as done in (|6ip . results 
in the following expressions for the component functions 



6 = 



(i-m+2) 



r 2 n' + Ar 2 WIi + 1671-/3(39 + p) 



+2iJ2 ( ~ e) VAU J 



Li 



A = - 



(J-l)G + 2) 



r 2 II + 4r 2 [/II 



2 % Y J ( - e) V A n J 



(70) 



(71) 



id 



One further differentiation of (|69[) gives an evolution 
equation for the axial sector of the Einstein equations, 
which at second order is 



1 

2^ 



(r 4 H) 



+ 2) 



n 



= 87re AB t, A \B+ie AB Y, ( ~ e)S A\B. 



(72) 



By expanding the energy-momentum perturbations in 
terms of the fundamental fluid variables (j46|). this equa- 
tion becomes 



2r A 



(I -!)(/ + 2) 



n : 



8/r 



8n(p + p)(3> - — [pCs' + v{p + p)]0 



AB 



(73) 



i.i 



The expanded form of this source is too large to be given 
here. 

In summary, the evolution of the second-order axial 
perturbations for Z > 2 is given by two coupled equa- 
tions; the matter equation (1641) determines the evolution 
of f3 and the master equation ((73)) that of the Gerlach 
and Sengupta scalar II. Reconstruction of the metric is 
facilitated by using Eqs. (|70H71|) . 



2. The case I = 1 

For I = 1 the right-hand side of the matter equation 
(|64[) slightly simplifies but the overall structure remains 



that of a transport equation for j3, 

(p+p)l$-cl(p+2U)P]= t J2{ { - £) l+^(r 2 ^£ A ) lA }, 

T,l 

,(74) 

whose sources are given by (|65|) and (|66j) particularizing 
to I = 1. However, the treatment of the gravitational part 
is rather different. The metric perturbation equation (|69[) 
is still valid for I = 1 and now becomes an equation for 

n, 



^e AB (r 4 n)l s = 87r^ + z^(- £ )^. (75) 



2r 



This allows direct integration of II from the knowledge of 
j3 and the sources, reflecting the non-existence of dipole 
gravitational waves. This is clearer by projecting the 
equation onto the frame vectors u and n A and expand- 
ing the fluid term according to Eq. (|46|) : 



-(n' + 4WOI) = -8n{p + p)f3-iJ2 i ~ e) 'PAU A 1 {76) 



1,1 



-(ri + 4[/TI) = -i^-^VAt 



(77) 



i,i 



In practice we plan to first integrate /3 using the evolution 
equation ([71)1 and then obtain II from the constraint (|76l) . 
The former implies dividing by the factor (p + p) which, 
at least for some equations of state, might vanish at the 
surface of a fluid ball but the solution can be shown to 
be nonsingular by power series expansions around the 
surface. 

A further difference is that now we cannot reconstruct 
ka from II using the perturbed Einstein equations. We 
need to integrate the definition (|6"51) of II. We proceed as 
follows. First, we note that any two-dimensional vector 
can be expressed in terms of two scalar functions y(x A ) 
and z(x A ) by writing 

~2^A = V\ A + £ A B z\ B = ~(y + z')u A + (y' + z)n A . (78) 

We insert this expression into the definition of the master 
scalar (1551) and obtain 



n 



-Z ]A \A 



Z + (AZ — VZ 



(79) 



Finally, we note that in the special case 1 = 1, the sym- 
metric trace free tensor X" 1 ^ defined in Eq. ((5|) vanishes, 
so that the condition M /i™ = in Eq. ([3]) no longer pro- 
vides a gauge condition. In RW gauge for / = 1 there thus 
remains one axial degree of freedom. The function y rep- 
resents this freedom and a gauge transformation y — > y+£ 
causes a change in the vector ka 



ka 



(80) 



Here the dots indicate source terms quadratic in first- 
order perturbations which are given explicitly in Ref. ■ 
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The resulting change in the component functions 6 and 
A defined in Eqs. ([701) and HHJ is 



6—>5- r 2 C + . 
A — > A + r 2 £ + 



(81) 
(82) 



We are therefore able to fix the gauge by eliminating 
either S or A or any combination of the two. It turns out 
to be most convenient to set A = 0. Equation (|78| then 
gives y' = — i, which allows us to integrate y from z and 
the values of y at r — (a residual gauge corresponding 
to a generator £ which obeys £' = 0). In this gauge, the 
vector ka is given by, 



ka 



-r 2 {y + z')u A , 



(83) 



where z is obtained from equation (|79|) and y' = —z. 



B. Polar perturbations 

Following decomposition (p32"j) . we can schematically 
write the second-order Einstein equations as follows. We 
indicate in each case the range of values of I, the label 
of the second-order perturbation, for which the various 
equations hold. For I > 0, 



Eab[IC] = -87r(/? + _p)7 + 47r( ( o - p)ip 



A B 

ABU n , 



,A„,B 



Eab [fc] = 8npuj + 8irp(rj — 4>) 

+ ^^ABU^, 

T,i 

n A n B E AB [K] = 8np{c 2 s uj + Ca) + 8tt P {t] + , 
+ J2 ie) VA B n A n B , 

T,i 



(84) 



(85) 



(86) 



E[K) = 8%p{cF s u; + Ca)+8TrpfC 

+ Y, {t) v- ( 87 ) 



For I > 1, 



l a E a [K] = -(p + p)a + J2 {e) T ? AU A , 

i,i 

i A E A [K] = J2 ie) VA nA - (89) 



i,i 



And finally, for I > 2, 



(90) 



The -E are the Gerlach and Sengupta linear differential 
operators which act on the polar metric perturbations 
{K.ab_J£} and are defined explicitly in Eqs. (D1-D4) of 
Ref. Q. 



1. The case I > 2 

By expanding the part of equation (|90[) which is linear 
in the second-order perturbations we arrive at the simple 
expression 



(91) 



1,1 



The scalar r\ is thus given directly in terms of the first- 
order perturbations. Note that, as commented at the 
beginning of this section, in the case of first-order per- 
turbation theory there are no source terms and the first- 
order r] vanishes. Even though, Eq. (|9"Tj) does not exist 
for I = 0, 1 and hence, in these cases, the first-order r] is 
nonvanishing. 

Motivated by the linearized equations of energy- 
momentum conservation, we can take linear combina- 
tions of the system (|84ti89|) such that it can be rewritten 

as 



I > : 

87r(p + p) 7 = (/C)' + C 7 + ]r«C 7 , 



1,1 

8irpu) = -K" + 2Uyj' + C u + (e) ^> 

T,l 

1>1: 

16Tr{p + p)a = V' + C Q +^ (e) C Q , 

i,i 

ii 



-x + x" + 2(/i - uW = s x + J2 {e) s x , 

T,i 

-ft + c 2 K" - 2c 2 s U^' = + (£) ^, 



(92) 
(93) 

(94) 
(95) 
(96) 
(97) 



i,i 



where C 7 , C u , C Q , S x , Sk: and are linear func- 
tions of the quantities Xt & and 7], their derivatives 
{x> x'> ^-j ^'j Vi v') v"} and of ifi and a. These sources 
are identical to those given in Ref. [45[ and are included 
in App.lAlfor completeness. We have also introduced the 
following linear combinations of quadratic sources, 



( £ »c 7 

(e)r? 

(e)C 
^ a 



^Vabu a u b 



"Pabu a u b , 



( e) S y = 2 (e) P + 4 



(0 



V A n J 



2 ie) V AB n A n B 



A{2v- W) {t) V A n A , 



^S K 



-c 2 u A u B +n A n B V< 



Vai 



4W (e) V A r A 
i n' 



-2^P A r^ 



(99) 
(100) 

3 

(101) 

(102) 
(103) 
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In order to illustrate the kind of quadratic terms that plicitly those corresponding to the last equation for the 
appear in the above equations, here we also present ex- variable tp, 



~ E 2 



' lT i^i{-^ + SX + fjtSX - fi5X + 55' + 55' - 2AA' - 2W55 + 2WXx\ 



+ Ef^A 2/C/C' + 2KK' + 2^(x + K)- 2(i/> + 2(m/j) ( X + K) - 167ra(p + p) (27 - ft) 



l(l + 1) 



3<5A -5X + 25X + pSX + 3p,5X + vXX + 55' + 355' - 4U5X - 2{3W - v)55 + (2W + u)XX 



+ 2( x + K) (i> + x' + K! + 2ni' - 2W X ~ 2WK) - 2ip(x + 1 + ft - 2Wft) \, 



(104) 



2i 
r 



2 E 2 _ x %\ X X + KX + X X + KX - # - H 



2iE\ 



i\tl\x{x + ^) - X + S (x + K) + 5£ + Xfo' + x") + n'ty - V (a + 5") + i/[3A(x + /C) - U 

2M ~- 1(1 + 1) -* l(l + l) - _ * _ - - - _\ 

— 5i> + y ^ 1 KX + v 2 2 y (O - xA + xA + 3/CA - + <fy>J - 8?rp^ 



5(x + K) ~2X4> 

+ 8tt( p + p) (27 -ft)f3- 4ttA(5p + 3p) (x + + 16ttA>(c^ + Cct) + £' 2^A - (/i - 2?7)l 
+ 6' + 1* ~ 2U(x + ft)] + (2U - /i)A> + 2(17 + p)Xft + A( M 2 + 3v 2 + 2U 2 ) (x + ft) 
+ (3v + 2U)j(X + S ft' -2(v-W)(x + ft) + (3u + 2W)Xx' + (3/xA - 5' + SUX + 2W5)k 
+ 2W 2 5rp + 2ipx(2v + W) - A(^x + t^ft + 2Wft) + 2v 2UXijj + 3WX(x + ft)- W5^ 
- 4 ftX- r 2 UW(6(x + ft) + Xft) + 2[i Udip + U(x + ft)X + 2(y + W)^X 

I 



(105) 



where, as before, we have also excluded the 1,1 = 0, 1 its expression obtained algebraically from the constraint 
cases by assuming vanishing fj and f] in order to avoid a 
larger expression. 

On the other hand, we note that the entropy perturba- f four evolution equations for the four unconstrained 
tion a is the only matter variable appearing in the linear variables x, ip and a, for I > 2. The additional con- 
sources for equations Q92J|97J). Its evolution follows from straint equations (|92II94|) provide algebraic relations to 



Therefore, Eqs. (|95H97I) and (|106[) form a closed system 



entropy conservation s = 0, 



<7+(7+fy = E (£)s - 



1,1 



reconstruct the matter perturbations 7, w and a. The 
energy-momentum conservation equations provide (re- 
(106) dundant) evolution equations for these other fluid vari- 
ables, namely, 



valid for I > 0, with the source terms 

- ~E^ l %{{2pX + 5X)s' -2aa) 
+ E° %{2j(iC + xW-(2j + ty 
+ (2fj-ft-x)d-}, (107) 
{ - ] S, = ~E^%{aXs'-0 + 5)a}. (108) 

The second-order variable 7 in (|106[) must be replaced by 



I > 

= + (109) 



1.1 



l + ^(7-|)-+cV = 5 7 +^^ 7 , (110) 



I > 1 

S a +Y / (t) Sa, (HI) 



1,1 
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where, as before, we have gathered the non-principal part 
in linear functions S u , S 7 and S a which are given in 
App. El We have further defined the quadratic sources 



P x 



w 5u = l l+Hue A u A + 2uMe 



r 2 ^8 AB 



,A _ (e) 



l A u A ,(U2) 



p x 



(p + p) x ^S a 



l l + V(*)E A n A + 2WV£ 



t 2 ^£ab 



n A - ^l A n A ,(U3) 



MS (i-m+2) (e)s 

2r 2 



1 



r 2(e) £ A + (e) X 
1^4 



(114) 



It is possible to use equations (|94|) and f|9T[) to eliminate 
derivatives of -0 in Eqs. (I109P and (|110j) . arriving at the 
sound wave equation 



1 + ^)7'=^ 



I > 1 



i,i 



i + y 7+ c y=^ 7 +E (e) ^ 



(115) 
(116) 



gauge conditions (H A = 0). In consequence, the vari- 
ables IC A b no longer serve as gauge invariants but in- 
stead change under gauge transformations generated by 
the vector £ M dx M = r 2 ^Z a dx a as 



K-ab ICab - {r 2 £\A)\B ~ (r 2 £\B)\A + 



K -> K - 2£ - 2r 2 v A £ l 



. (H7) 
, (118) 



where the dots indicate source terms quadratic in the 
first order perturbations. 

At this point, we follow Campolattaro and Thorne [34| 
and fix the remaining gauge degree of freedom by de- 
manding K, to vanish. Even then we are still left with the 
residual freedom of transformations under the restricted 
class of functions £ such that £ + r 2 v A £i A = 0. 



We have already mentioned that Eq. (I91j) is absent for 
1 = 1, and hence we proceed differently than for I > 2. 
We use equations ([TPS]) . (ITTTj) . ([TT51) and (fTTB|> to evolve 
the matter perturbations {a, a, oj, 7}. These equations 
contain the metric perturbations {T],tf),x} but n ot their 
derivatives. On the other hand, we obtain the metric 
perturbations from the Einstein equations (|92M94I) and 
(|96H97[) . With JC = from the additional gauge free- 
dom, these five equations contain at most first deriva- 
tives of {rj, ip,x}- Hence, there are five equations for six 
unknowns: we can obtain {t/j, ip',x,x'} but only the fol- 
lowing combination of r\ and rf , 



Again, the linear sources S u and 5 7 are given in App.lAl 



Dr] 



V\A 



1 



r\ B r\ B 



■(yVrf - Urj), 



(119) 



2. Polar perturbations with I = 1 

In the special case I = 1 there are still three polar de- 
grees of freedom, but the RW gauge only imposes two 



where |i>| 2 = v v A . In consequence, we can only inte- 
grate j] on a spatial surface that is everywhere normal to 
r = const, surfaces, and therefore we also integrate the 
other scalars tp and x m the same form, 



r\v\ 2 D V = rhs(Al2) - i £ {^S K + (1 - c 2 s )^C u - 2U^C a }, (120) 

T,i 

r\v\ 2 D X = rhs(Al3) + ^{2[/( £ )C Q - (£) C ti; }, (121) 

1,1 

r\v\ 2 D^ = rhs(Al4)-^{( e )C7 7 + W (e) C a -C/^^}, (122) 



where the notation rhs(X) stands for the right-hand side 
of Eq. (X) of Ref. j4j|. 



3. The case I = 



The / = perturbations are treated in a manner simi- 
lar to the case 1 = 1. We evolve the matter perturbations 
and then obtain the metric perturbations from the con- 
straints. This is a direct consequence of the absence of 
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radiative degrees of freedom in spherical symmetry. 

For I = we still have four perturbative variables, 
namely x, 77, -0 and JC, and only two gauge degrees of free- 
dom in the gauge-generator vector. The RW gauge does 
not impose any condition for / = 0, and hence we need to 
fix those two gauge degrees in some other way. We do this 
by extending the gauge used in the previous subsection 

* = °> * = pn^fo-x)- ( 123 ) 

In polar-radial coordinates the second condition implies 



J 



the vanishing of the variable ip. A further feature of this 
case is that the velocity perturbation a also vanishes. 

The evolution procedure is then as follows. The mat- 
ter perturbations a, uj and 7 are evolved by Eqs. (|106|) . 
(|109|) and ()110|) . respectively. The latter two equations 
contain derivatives of the metric perturbations that can 
be removed by using the perturbed Einstein equations 
(|84II86|) . The resulting equations, as well as the constraint 
equations for the two non-vanishing metric perturbations 
rj and \ read as follows, 



(1 ^J-' 



rhs(A15) + £ {&S U - { Jt±A{Uu A u B - Wu A n B )^V AB }, 



l,l 



l + -)7 + cV = rhs(A16) + ^{ (e) S 7 



l,l 



2p\ 

(P + P) 
2p\v\ 2 



{Uu A n B -Wn A n B )^V AB }, 



r\v\ 2 Dr] = rhs(A18) + -^V {{U 2 + W 2 ){u A u B + n A n B ) - 4UWu A n B \^V ABl 
2mr — ' 



l,l 



r\v\ 2 D X = rhs(A19) + t^^U^ 2 + W 2 )u A u B - 2UWu A n B }^V ABl 



(124) 
(125) 
(126) 
(127) 



l,l 



where we have again used the notation rhs(X) to refer to field normal to the surface is defined as 

?V = <pf,n, where ip = (/,„/'") 



the right-hand side of different formulas of Ref . [45 1 



-1/2 



(128) 



V. PERTURBATIVE MATCHING TO AN 
EXTERIOR VACUUM 

So far, we have studied the interior of a perfect fluid 
system. We now assume that our system is a fluid 
star surrounded by vacuum, with both regions, interior 
and exterior, separated by a surface £ where the pres- 
sure vanishes. Hence, the background exterior will be 
Schwarzschild. 

The first part of this section describes high-order per- 
turbative matching across any timelike surface in a gen- 
eral background spacetime. The second part will partic- 
ularize to the case of a spherical fluid interior matched 
to a vacuum exterior. This generalizes the first-order re- 
sults of Ref. [HI for the same scenario, and we closely 
follow their notation. 



A. High-order matching conditions 

We describe the matching surface S as the zero level 
set of a smooth scalar field f{x >l ), with arbitrary smooth 
continuation off the surface. The unit spacelike vector 



[lis — 

i a „ of 



From this vector we construct the induced metric i 
9n» — n^n v and the extrinsic curvature = 
the surface S, both of which must be continuous across 
S to ensure a smooth matching between the interior and 
exterior solutions. Before discussing high-order perturba- 
tions of these objects we need to address two important 
issues arising in the perturbative version of the matching 
problem: index positioning and gauge dependence. 
For a generic tensor field T M we have 



A[T M ] + g^A[T% 



(129) 



because in general the perturbation of the metric field 
does not vanish. Imposing continuity on the perturba- 
tions of the covariant or contravariant forms of the ten- 
sors i and e can therefore lead to different results, and we 
must decide in advance which are the proper quantities to 
be used. The first-order discussion in [54[ shows that we 
must use perturbations of the contravariant fundamental 
forms. Essentially, this is because a contravariant ten- 
sor field T^ v - - is intrinsic to the surface S if and only 
jf T^ v ' f jfl = on any of its indices /i. The equivalent 
condition for a covariant tensor field T^..., on the other 
hand, would be T^ u ,,,g^ a f_ a = 0, which involves the am- 
bient metric and therefore introduces additional informa- 
tion not intrinsic to the surface. This argument gener- 
alizes to higher-perturbative orders, as we will later see, 
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and so we will impose continuity of A™ [i^] and A" \e^ v \ 
for all n. 

Second, we need to deal with the gauge freedom aris- 
ing from the arbitrariness in our choice of mapping 
$ between the perturbed and the background space- 
times. Under a general perturbation the scalar / will also 
change, so that the perturbed surface will be described 
by the level surfaces of / + A$[/] + ... (a field on the 
background manifold), where we have made explicit the 
gauge $ relating the background and perturbed space- 
times. It is possible, therefore, for a perturbed interior 
point to be located at a coordinate position correspond- 
ing to the background exterior. Such a situation can 
be handled consistently using one-sided derivatives f5Qj j . 
There remains, however, the question of the correct con- 
tinuity conditions, because these conditions may not be 
equivalent if expressed in different gauges. A convenient 
treatment of this difficulty is possible because there exists 
a privileged class of gauges fi, characterized by the con- 
dition A q \f] = which generalizes the first-order surface 
gauge of [54[. This does not imply that the shape of the 
surface will not change; indeed, the surface can be highly 
distorted in the perturbed manifold. However, this kind 
of mappings between perturbed and background space- 
times will identify any point of the perturbed surface 
with another point of the background surface. In surface 
gauge, the matching conditions at any perturbative or- 
der are therefore given by the continuity of the induced 
metric A^[z M "] and the extrinsic curvature A^[e F ]. 

We emphasize that surface gauge is only used to de- 
fine the continuity conditions. We can still work using 
any other gauge as convenient for the interior or exterior 
problems. The gauge freedom is then handled by con- 
structing gauge- invariants associated with surface gauge, 
that is, combinations of the perturbations in an arbitrary 
gauge whose values coincide with the result in surface 
gauge. See [2] for a discussion on the construction of 
such gauge-invariants. This can be achieved by finding 
the general form of a gauge transformation from arbitrary 
to surface gauge. Such transformation will be parame- 
terized by the gauge vectors { fl} £, defined by 
solving the following equation order by order, 



= A"[/]+E 



1 



- (n-m)! ( |- ) 2!^7.(m!) fc ™fci!...fc m ! 

xC% r X k { Z }i A n - m [f}. (130) 

This equation is just the gauge transformation (|TU)) of 
the n th order perturbation of the scalar field / to surface 
gauge. In particular, solutions at first and second orders 
are given by 

<% = -A[/K, (131) 
<% = -A 2 [/]rv + /A[/]A[/], M . (132) 

These solutions fix only one of the four degrees of freedom 
in each generator. There thus remain three additional 



degrees of freedom which represent gauge changes among 
different surface gauges. 

Now we can define the gauge-invariant combinations 
of perturbations via an operator A whose action on a 
background tensor field T M is defined as 



S"[T"] = A n pT"] +^ -7— 

771=1 ^ 



E 

n) 

A"- m [T^],(133) 



to)! ' fci \...km ! 

(K m ) 



2! fe 2...( m !)fc™ {1 >?' {m >« 

where the gauge vectors { are those obtained 

by solving equation (|130[) . Note that the gauge genera- 
tors depend implicitly on the metric perturbations, so 
that this formula is non-linear and highly nontrivial to 
implement. 

The A operator has been constructed explicitly so that 
it obeys A / = for all n. Specifically, this implies 



A' 



5>1 



For example, at first order we obtain 

A>] = |ri Q [n* {li h a0 - 2(^A[/]), a ] 



(134) 



(135) 



We do not display explicit second- and higher-order for- 
mulas here as they are rather complicated and do not 
contribute to enlighten the discussion. 

Equation (|134j) implies that perturbations of any con- 
travariant tensor intrinsic to the surface E will also be 
intrinsic to the perturbed surface. In other words, for a 
background tensor with T^n^ = 0we also have 



A 



Q. 



(136) 



Note that covariant tensor fields do not share this prop- 
erty. 

In summary, we require the barred perturbations of the 
contravariant fundamental forms A [i^ v ] and A [e^] to 
be continuous across the surface E at any perturbative 
order n. Note that these conditions are formulated in 
any gauge. In contrast to Ref. [HH we therefore do not 
need to impose surface gauge. 

We conclude this discussion by mentioning the alter- 
native approach to perturbative matching introduced by 
Mukohyama [6(|] and further developed in [6lJ . Their re- 
sults have been given for first-order perturbations and co- 
incide with those presented here and in [54[ . Mukohyama 
and coworkers do not impose surface gauge and include 
the gauge freedom within the matching surface E. Com- 
pared with our approach, the main difference is their use 
of an abstract copy of the surface E. The matching is per- 
formed separately between the boundaries of the interior 
and exterior spacetimes and that new surface. This pro- 
cedure introduces an additional geometric structure and 
gives rise to a new type of gauge invariance and, thus, 
the concept of double gauge-invariants. While this can 
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be a convenient feature, as for example in the reduction 
from a D-dimcnsional spacetime to a (D — l)-dimensional 
brane, it would complicate our comparatively simple sit- 
uation of a spherical background, where the geometry of 
the matching surface is trivial. 



B. Matching to vacuum 

In spherical symmetry, the matching conditions can be 
decomposed into tensor spherical harmonics, which al- 
lows us to extract or inject information through the sur- 
face independently for the different harmonic components 
of the metric and matter fields. We have constructed 
those decomposed continuity conditions for second-order 
perturbations, but we will not discuss them here, be- 
cause they involve very large expressions. Again, they 
coincide with the expressions given in adding sources 
quadratic in first-order perturbations. In App. |B]we de- 
scribe in detail and explicitly how to perform the match- 
ing for the particular case of a first-order (Z = l,m = 0) 
axial and a second-order (Z = 2, m = 0) polar mode. Here 
we will also mention different ways of describing the ex- 
terior perturbations. 

The natural exterior frame is given in terms of the 
radial vector and its orthogonal defined as 



.1-4 



-e AB r B , 



(137) 



where the normalization factor a is given in terms of the 

Hawking mass as a -1 = yjl — 2M, l n the interior of the 

star these vectors are related to the fluid frame (u A ,n A ) 
by a hyperbolic rotation. 

The axial matching is simplified by the existence of the 
Gerlach and Sengupta master scalar, which can be de- 
fined both in the interior and the exterior withoug using 
fluid information, and so it obeys very simple continuity 
conditions. The polar problem is harder because none of 
the internal perturbations matches easily with the natu- 
ral variable describing vacuum perturbations, the Zerilli 
scalar 



2a- 2 r 2 



6M + (Z-l)(Z + 2)', 



■(r B ™K AB -ar™K\ A 



(138) 



This function obeys the following unconstrained evolu- 
tion equation [62j], 



{ ">*| A A -U M * 



(139) 



where V is a potential that depends on the scalars (r, M) 
and on the harmonic label Z, whereas {ni S is a source 
that depends on lower-order perturbations and vanishes 
for linear perturbations n = 1. At second and higher- 
orders one must deal with unexpected divergencies at 
asymptotic spatial infinity, which must be regularized by 
modifying the definition of the Zerilli scalar with ade- 
quate lower-order terms, and hence also modifying the 



form of the source of equation (|139[) . See [62j for details 
on the regularization procedure. 

Using the gauge invariance of our perturbative func- 
tions, we could use a particular gauge to simplify some 
equations. For instance, imposing a gauge given by the 
conditions G = 0, K = 0, and r A r B H AB = 0, the Zerilli 
scalar is essentially the harmonic coefficient r A H A , 



BCL 



2er 1 H(Z + l) 
~6M+(Z-l)(Z + 2)r 



(140) 

with lni R being a nonlinear source depending on lower 
order perturbations. The mentioned gauge is similar to 
the one introduced in Ref. [63] (G = 0, K = 0, and 
t A t B H A b =0). In fact, making use of the perturbative 
field equations, the linear part of the Zerilli function turns 
out have the same form. 



VI. CONCLUSIONS 

In this work we have presented a general formal- 
ism for second-order perturbations of a time-dependent, 
self-gravitating, spherically-symmetric perfect fluid. To 
this end we have extended the formalism developed in 
Refs. [l|, [2[, which in turn can be viewed as a gener- 
alization of the Gerlach-Sengupta formalism [43|, |44| to 
higher orders. The starting point is a 2+2 decompo- 
sition of the background spherical spacetime; the four- 
dimensional manifold M. is given by the product of a 
two-dimensional Lorentzian manifold M 2 and the two- 
sphere. This separation enables us to fix the standard 
angular coordinates (6, <f>) on the sphere while formulat- 
ing the fields on 2 in a manifestly covariant form. Af- 
ter decomposition into tensor spherical harmonics we are 
left with tensors that depend only on the coordinates 
on M 2 . We have further defined gauge-invariant objects 
(for Z > 2) in a manner which identifies them with the 
quantities obtained in RW gauge. 

An important feature of the perfect fluid in spheri- 
cal symmetry is that it defines a complete basis of vec- 
tors on the manifold Ai 2 , given by the four- velocity of 
the fluid and its orthogonal. This basis has been used 
to project all tensorial, background and perturbative, 
quantities in such a way that one works exclusively with 
scalar fields. For the case of second-order perturbations, 
we have explicitly decomposed the perturbations of the 
energy-momentum into tensor harmonics. By doing so, 
we have presented for the first time the complete matter 
sources that appear in the perturbative evolution equa- 
tions. 

Next, we have combined the geometric sources for the 
Einstein and matter equations, as given in Ref. [l|, with 
the matter sources derived in this article and, thus, ob- 
tained the complete set of evolution equations. This has 
been achieved by following the linear analysis of Ref. [45[ 
for both the polar and axial modes and arbitrary har- 
monic labels Z. Except for the polar case with Z < 2, we 
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have thus obtained a combined set of evolution equations 
and constraints. The metric perturbations with I < 2 
have no radiative freedom and are governed exclusively 
by constraints. In this solution is obtained by first 

evolving the matter perturbations and then obtaining the 
metric perturbations from the constraints. 

We have also analyzed the problem of matching the in- 
terior perturbative perfect fluid solution with the exterior 
vacuum spacetime at any order. We have demonstrated 
why tensors are naturally matched in their contravariant 
form. In surface gauge, defined as the gauge where the 
perturbed surface is mapped to that of the background, 
the perturbations of the contravariant fundamental forms 
must be continuous across the surface. Even though this 
result is valid as such in surface gauge only, the relations 
can be transformed to arbitrary gauge so that practical 
applications are not necessarily restricted to the use of 
surface gauge. 

In summary, we have presented a general formalism 
capable of studying coupling between any two first-order 
modes of a spherical time-dependent star. In future work, 
we plan to apply this formalism to a variety of stellar 
models including, but not restricted to, the emission of 
gravitational waves by slowly rotating, collapsing stars. 
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Appendix A: Polar linear sources 



In this appendix we give the linear sources appearing 
in the polar sector of the evolution equations (|92l) - (|116l) . 
These sources have been already given in Ref. |45|, but 
have been recomputed in the course of this investigation, 
as an intermediate step towards the construction of the 
second-order sources. There is full agreement with the 
results of Ref. El. 



2i> 2 + &wp-^-- 2U(p, - U) 



(!-!)(/ + 2). 



3px + 4(a* - U)t - (5u - 2W)x' - 2[2p,v - 2(p - U)W + p! - z>]V> 
2rj" - 2(p - U)r) + (8v - 6W)r]' - 



Au 2 + '( f + 1 ) + 8 + 8vW + A(2 P U + U 2 - AW 2 - 8tt P ) 



(1 + c 2 s )U x + W + c 2 (p + 2U)]K - W(l - c 2 s ) x ' -{v + 2W(? S )K' 



r - W 2 
2r 2 

2Uf] + 2Wr]' 



8np - 



2 fKi + i) 



+ 2U(2n + U) - 8tt P 



(x + ic) 



(1 + c 2 s )x + 2[-p,W(l - c 2 ) + {v + W)U(l + c 2 s )]i> + 8irCpa 



1(1 + 1) + 2 



6W 2 + 167rp - 2U{2 P + U)c 2 s 



V, 



2v{ x + K) + 2/iV> + x' ~ 277(1/ -W)- 2r)', 



C 7 = -Wx + Ux- {p-2U)JC' + 



1(1 + 1) + 2 



T), (Al) 



1(1 + 1) 

v 2 ; + 2U(2p, + U) - 8tt P 



+ 2U(2 P + U)- 2W(2v + W) + 8ir(p - p) 



U + K)- {l 1) { [ + 2) X + 2[uU + + U)W]i> 



U X +(^ + 2U)t + Wx - 2WK' - 2r]U(2p + U), 
2/i(x + K) + 2vi) + x + 2K - 2r)(fi + U), 



(A2) 
(A3) 

ijj - 2Ur)', (A4) 



(A5) 
(A6) 
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Suj — 



1(1 + 1) x + 3K. 
n — a H 



U + 2W- 



7+77 



+ 0* + c 



C 



See — 

Sui = 



i + £ 
p 

Ca' -a 
1C + X 



2 / ,; 

x' + £' - 2,?' 



+ [ c 2( At +2C7)- M ] [j-t 



, 2 P P 

v id 



pdc 2 s 



P 



dp 



c? ds 



,dC 



+ i] - c? s (n + 2U)a + 



ds 

duj + Ca 



P/ c s ^ 



M_cfl + — — 
ds \ <? s dp 



1 



(A7) 

(A8) 
(A9) 



. l Jl±R + 8ir(p + p) ] a + £ + (m + C/)^ - M ( x + £) 



+ (/* + 2J7) c 



w 



C*(/j + 2t/)o- 



i/ 1 



(A10) 



P 



KJ 



+ (cl(p + 2U)- fJ ) 7 



c 



p 



+ s' 



/i^ - v(x + K) + {v - W)r) 



1M 

c? ds 



, P 2 

^ c„ 



(All) 
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Appendix B: Example of matching: first-order axial 
(1=1, m=0) and second-order polar (1=2, m=0) 
modes 



The general expressions for the matching at the stellar 
surface are quite large. In order to illustrate the pro- 
cedure in detail, we therefore analyze in this appendix 
the matching conditions for the particular case of a first- 
order {/ = l,7n = 0} axial mode that, by self coupling, 
generates the second-order {I — 2,m — 0} polar mode. 
This is a particularly interesting situation since the non- 
radiative first-order mode can be interpreted as a slow 
rotation of the star that produces gravitational radiation 
through self coupling. For identification of the pertur- 
bative order of the individual terms, we use a left su- 
perindex {1} on first-order perturbations. 

The background junction conditions are straightfor- 
wardly deduced from continuity of the two first funda- 
mental forms of the surface. The continuous quantities 
include the scalars r, v, and W defined in (f2"Tj) and (j2li|) 
respectively. Derivatives of continuous quantities in the 
direction of the fluid velocity u A must also be continuous. 
This leads to continuity of U (l23|) and, since both deriva- 
tives of r are continuous, to continuity of the Hawking 
mass (gSJ). 

In our particular case, the negative pressure of the fluid 
—p will be interpreted as the scalar function /, since the 
stellar surface is characterized by p = 0. The negative 
sign has been chosen such that / increases with radius r. 
Furthermore the pressure must be continuous across the 



surface, whereas the energy density p may jump there. 

The first-order axial matching is simplified by the ex- 
istence of the Gerlach and Sengupta master scalar ([551) . 
which can be defined both in the interior and the exte- 
rior without using fluid information. In fact, a first-order 
junction condition for our particular example is directly 
obtained from the continuity of the master scalar {1 TI 
across the surface. Regarding the axial vector {1} ka, we 
only have continuity of the timelike component {1} kau a , 
that has been defined as {1} 5 in Eq. (|6Tj) . In the gauge 
proposed at the end of Subsec. HV A2| {1} A = {1} n A n A = 
0, however, the entire axial vector turns out to 

be continuous. For simplicity and consistency with that 
subsection, we will employ this gauge here. A further 
quantity that must be continuous is given by the combi- 
nation 



^U' + ^p^p 



(Bl) 



which depends on the axial fluid perturbation {1} /3. 

The second-order polar problem is more complicated 
because none of the internal perturbations matches 
straightforwardly with the natural variable describing 
vacuum perturbations, the Zerilli scalar. Decomposing 

2 2 

into harmonics the objects A [i^] and A [e MI/ ], we find 
that the following second-order polar quantities must be 
continuous, 



At = N + S, 

A 2 = X + 2(i/ + W)N - 2?7, 



(B2) 
(B3) 



19 



A 3 
A A 



K - 2NW, 

ij) + 2N - 2UN, 

K! - 2t]W + 2UN + 2 \4%p - 2U 2 + W 2 



vW - 



M + 4r 



V, 



(B4) 
(B5) 



(B6) 



A 6 = x' + 2fiijj-2r]' + 2(W + v)r)-2N + 2UN 

( M — 5r\ 
- 2 4ttp+ 2v 2 + U 2 -vW+ 5 — N,(B7) 



where the quadratic source 

s = 1 



5?r 



{1}^{1}^ 



(B8) 



appears only in the first continuous combination A\ , and 
•q is given algebraically in terms of first-order perturba- 
tions (Bp , 



r, = -^{%ir( P + p) {1} p 2 + \[{v-2W)^5+ (1} <5'] 
V57T I 2 



(B9) 



Note that some terms of the right-hand side of this ex- 
pression are continuous, e. g. {11 S 2 /r 2 . Hence, when 
introducing this expression in the definition of the con- 
tinuous objects (|B2IIB7[) . those terms can be removed. 
We note, however, that special care is required for the 
term rj' in Eq. (|B7[) , because prime derivatives of con- 
tinuous objects do not have to be continuous. With the 
exception of N, all perturbative objects that form part 
of the expressions for the continuous objects have been 
defined in the main body of the article. The remaining 
variable N is proportional to the gauge invariant associ- 
ated to the pressure perturbations p by 



<P 



■2£ {1}p A[ P }+ £ {2 



(BIO) 



The pressure p must not be confused with the vector 

whose harmonic components are given in Ref. 
Last relation can be written in terms of the second-order 
gauge-invariant perturbations of the energy-density and 
entropy, 



N 



-<pp{c 2 uj + Co). 



(Bll) 



The subindex GI denotes the perturbation expressed in 
terms of the gauge-invariant objects that, again, has a 
form equivalent to the perturbation in the RW gauge 
Mpfi _ q This last equation has no quadratic terms 
in first-order energy-density and entropy perturbations 
because both are polar and we assume the first-order 
perturbations to be exclusively axial. The norm of the 
normal vector is defined by (|128[) . 



p = {pap ) 



1 



(B12) 



The second equality holds because the pressure vanishes 
at the surface at all times; therefore p = 0. The minus 
sign comes from the fact that p 1 < 0. On the other hand, 
making use of the background Euler equation for the fluid 
(f30|) . we obtain 



P = -v(p + p). 



(B13) 



Combining the last three equations, we finally obtain N 
in terms of the fluid variables on the surface, 



N = - 



c 2 uj + Co 



(B14) 



As expected, the continuity conditions (|B2HB7[) reduce 
to the first-order expressions given in [54j in the absence 
of the source S and the component r). 



Extraction 



In the interior of the star, the relation between the 
natural exterior frame (r A ,t A ), defined in (|1 3T[) . and the 
fluid frame (u A ,n A ) is given by a hyperbolic rotation, 



if 



-arUu A + arWn A , (B15) 
-arWu A + arUn A , (B16) 



Replacing this form of the radial vector r A in the defini- 
tion of the Zcrilli function (| 138[) and rewriting it in terms 
of the fluid variables, one obtains 



2r 4 



-TT- l^out + 2UWi/j out - WJC' m 
1 r l 



6M + (I + 2) (I - 1) 5 
2U 2 7] out + (U 2 + W 2 ) (xout + /C out )} + r/C out , (B17) 



where primes and dots are always expressed in the fluid 
frame, and the expression is evaluated just outside the 
surface. Making use of the continuous quantities (|B2l 
IB7[) , we arrive at the formula that gives the outside Zerilli 
function in terms of the fluid inner variables, 



2r 4 



WICL 



6M + (l + m -l)A U,C * + 2UW ^ 
2U\ n + (U 2 + W 2 ) (xm + Kkn) - ^r 2 P WN; 

w 



-(-8 + / + / 2 )(S' out -5 in )}+r/C il 



(B18) 



At first order, because of the vanishing of the source 5, 
the last term in curly brackets would disappear. 



2. Injection 

Inside the star, the polar variables x ancl K satisfy a 
wave equation. Therefore, boundary conditions on the 
surface of the star must be given for these variables. 
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The tensor K-ab is decomposed in Schwarzschild coor- 
dinates as, 

^-(£: f:)- ( bi9 ) 

Outside the star the polar metric perturbations 
{H rr , H rt , H tt ,JC] can be recovered once the Zerilli func- 
tion has been determined (62[. Transfer of information 
from the exterior to the interior of the star is again ob- 
tained from a hyperbolic rotation and use of the continu- 
ous quantities (IB2HB7[) . Specifically, using the continuous 
quantities (|B3[) and (|B4[) . we get 

Xin = Xout - 2(l/ + W)(Sout - Sin) 



- 2(? hut ~ mn ), (B20) 
Kin = K out + 2W(S out -S in ). (B21) 

Applying a hyperbolic rotation results in 

Xout = (r _ \ M)2 {rWHtt + 2r 3 (r - 2M)UWH rt 
+ r 2 (r-2M) 2 W 2 H rr } -/C ou t- (B22) 

Hence, we have succeeded in writing the interior 
variables {IC,x} in terms of the exterior variables 
{IC, Htt, H tr , H rr }, which, in turn, can be obtained in 
terms of the Zerilli function. 
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